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SUNNARV 

A constitutive theory Is presented for representing the anisotropic v!$- 
coplastlc behavior of high-temperature alloys that possess directional proper- 
ties resulting from controlled grain growth or solidification. The theory Is 
•n extension of a viscoplastic model that has been applied In structural analy- 
ses Involving Isotropic metals. Anisotropy Is Introduced through the defini- 
tion of a vector field that Identifies a preferential (solidification) direc- 
tion at each material point. Following the development of a full multlaxlal 
theory, application Is made to homogeneously stressed elements In pure shear 
and to a uniaxial 1y stressed rectangular block In plane stress with the stress 
direction oriented at an arbitrary angle with the material direction. It Is 
shown that an additional material parameter Introduced to characterize the 
degree of anisotropy can be determined on the basis of simple creep tests. 


INTRODUCTION 

The need for greater efficiency In aircraft engines places Increasing 
demands on the high- temperature structural alloys used for engine components. 
As higher operating temperatures are sought, advanced materials are being de- 
veloped to meet these Increased demands. Good examples are the single crystal 
(SC) and directionally solidified (DS) polycrystal line materials finding ap- 
plication as turbine airfoil components. An advantage of these materials over 
conventionally cast alloys Is their Increased strength (e.o., creep and creep- 
rupture strength, yield strength, etc.) In a preferential (grain orowth) di- 
rection, which In the case of a turbine blade can be advantageously oriented 
radially (centrlfugally). Improved creep and creep-fatigue properties result 
as well as reduced susceptibility to grain boundary corrosion and oxidation. 

The directional properties of SC or DS metals render them highly aniso- 
tropic relative to conventional alloys. This Introduces additional complexity 
In understanding and mathematically representing their mechanical behavior 
over and above the already enormous complexities associated with elevated 
temperature. 

Here, the unified constitutive model of Robinson (refs. 1 and 2) that has 
found application In representing Important behavioral features of high- 
temperature Isotropic metals Is extended to account for the effects of aniso- 
tropy. Each material point Is taken to have a uniquely Identifiable direction 
designated by a vector. An extended material body Is thus treated as being 
locally transversely anisotropic although the preferential direction may vary 
from point to point as represented by a vector field. It Is believed that 
this relatively simple model captures the essence of anisotropy as Induced by 
directional grain growth and solidification without undue complication. 

‘University of Akron, Akron, Ohio and National Research Council - NASA 
Research Associate. 


The Isotropic, Isothemel for* of the Robinson vlscoplestlc theory Is 
first discussed with emphasis on Its derlvablllty from • potentlel function. 
Full Isotropy Is treeted by taking the applied and Internal stress dependence 
of the potential function In terms of the principal Invariants of the stress 
tensors. The extension to anisotropy Is made by replacing the principal In- 
variants with another set of stress Invariants that reflect the appropriate 
material symmetry. 

Following the development of the full multlaxlal theory, application is 
made to simple states of shear stress oriented transverse to and along the 
preferential material direction. A final application Is made to a uniaxial 1y 
stressed rectangular block of material In plane stress with a uniformly or1en< 
ted material direction taken at an arbitrary angle with the direction of 
stress. 


SYMBOLS 

a^ j components of devlatorlc Internal stress 

dj components of unit vector 

F scalar function of stress 

f material function 

6 scalar function of stress 

g material function 

g^ hardening function 

g r recovery function 

H material constant 

h hardening function 

1^ Invariants of effective stress 

invariants of Internal stress 
0^ principal Invariants of effective stress 
principal Invariants of internal stress 
K threshold transverse shear stress 

K d threshold longitudinal shear stress 

m material constant 

n material constant 

R material constant 

s Internal stress In transverse shear 

I internal stress In longitudinal shear 
components of applied devlatorlc stress 
x^ coordinate directions 

a uniaxial Internal stress 


2 



$ material constant 

t shear strain rate 

a.|j Kronecker delta 

c uniaxial Inelastic strain rate 

tjj components of Inelastic strain rate 

o potential function 

a uniaxial normal stress 

components of applied stress 
t shear stress 

Xjj components of effective stress 

f angle between Xj axis and stress direction 

u material constant 
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THE ISOTROPIC VISCOPLASTIC MODEL 

The flow and evolutionary equations In the Robinson model are taken to be 
derivable from a potential function q of the applied and Internal stress. 

The components of these stress tensors are denoted by ojj and ~ij, respec- 
tively. Thus, we have 

fi ■ j»°i j)« U) 


For the sake of simplicity, the present development Is restricted to iso- 
thermal conditions. Extension to nonisothermal conditions follows the devel- 
opment presented in references 1 and 2. 

As moderate hydrostatic stress is known to have essentially no effect on 
inelastic behavior, the stress dependence is taken in terms of the devlatoric 
components of the applied stress 

S 1j " °ij ‘ i °kk a 1j * 2 


and of the internal stress 


a 


ij 


®1j “ 1 “kk°1j 


1 


( 3 ) 


where the symbol j denotes the usual Kronecker delta. We further Identify 

“ S ij ~ a ij 

as the effective stress. The potential nature of o is expressed by 


c 1j 



( 5 ) 
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and 
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where 4. . represents the Inelastic strain rate and h Is a scalar function 
of the Internal stress. Equation (5) Is termed the flow law. Equation (6) Is 
termed the evolutionary law. 

The appropriateness of equations (5) and (6) has been discussed on physi- 
cal and thermodynamical grounds by several authors Including Rice (ref. 3), 
Ponter and Leckle (ref. 4), Valanls (ref. 5), and Robinson (ref. 6). Equations 
(5) and (6) are shown In reference 4 to hold exactly for an Individual slip 
system In a polycrystal line metal deforming at high temperature where Is 
Interpreted as the local Internal flow stress on the slip plane. Equation (5) 
remains valid (ref. 3) for a polycrystalllne metal where o.. and 4., are 
Interpreted as the average stress and Inelastic strain rate ‘over a volume of 
material that Is large compared to the crystal size. The derivative In equa- 
tion (6), however, as applied to an Individual slip system, requires constant 
local stress, whereas In the present context It Implies constant average 
stress. As pointed out In reference 4, It Is not generally possible to assume 
that constant average stress Implies constant local stress and, consequently, 
that equation (6) remains exactly true for a polycrystalllne material. Nev- 
ertheless, constitutive relationships have been derived from equations (5) 
and (6) that are consistent generalizations of well established classical 
equations (refs. 3 and 4) and that are capable of accurately representing 
Important features of high temperature behavior of metals Including rate- 
dependent plasticity, creep, recovery and their Interactions (refs. 1 and 2). 

The function a In the Robinson model can be written as 

a - K 2 (y^^(F)dF ♦ J"\ g(G)<H>l (7) 

where the stress dependence enters through the scalar functions 

F(r i j) and 6(a iJ ) (8) 


taken as depending on the effective and Internal stress, respectively. The 
functions f and g and the material parameters K, u, R, and H are 
assumed known for present purposes; they are determined as described In 
earlier writings (refs. 1 and 2). 

Under full Isotropy, the functions F and G can be taken to depend on 
the principal stress Invariants 


and 


1 > 
°2 " 1 £ 1j E ji 

► 

J 3 “ 1 I 1j I jk £ k1 


(9) 
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SfftSW 


r 2 m l a 1j a J1 


h “ T a u a jk a ki 


( 10 ) 


In the spirit of von Mises, we retain only J2 and fi% Quadratic In 
stress , and take 


F ■ -y - 1 

ir 


(U) 


and 




(12) 


Equation (11) plays the role of a (Bingham) yield condition with K de- 
noting the magnitude of the threshold shear stress below which Inelastic defor- 
mation does not occur - inelastic strain occurs only for F > 0. For our pur- 
poses, we treat K as a constant; more generally It is considered a scalar 
state variable . 

1 he flow and evolutionary equations are determined directly from equations 
(5) and (6) making use of equations (7) to (12) and taking 


H 

gJET 


(13) 


in equation (6). The details of this development are given In appendix A. Here 
we state just the result, i.e., the flow law 


2« f j - f(F) tlj 


(14) 


and the evolutionary law 


*1J “ gJGT *ij ' Rg r (G > a ij 


(15) 


in which 


• 9 r (G) - (16) 

This is essentially the form of the Robinson model for a fully Isotropic 
material and for isothermal conditions. Some Important features of the model, 
such as the accompanying inequalities (refs. 1 and 2), are not expressed or 
discussed here as they do not pertain directly to the extension to anisotropy. 


5 



The evolutionary law (eq. (15)) is of the widely accepted Bailey-Orowan 
type, which presumes that high- temperature deformation occurs under the action 
of two simultaneously competing mechanisms, a hardening process proceeding with 
accumulated deformation (characterized by the first term in eq. (15)) and a re- 
covery term proceeding with time (characterized by the second term). Steady 
state then corresponds to the condition where the two competing mechanisms 
balance and h. . » 0. 

' V 

In most applications of the theory to date, the function f has been 
chosen as 


or 


cr:g;^a ? . pact \s 
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f(F) - F n 

► 


(17) 


f (F) « (sinh F) n J 


and, as suggested by the experimental results of Mitra and McLean (ref. 7), 


g(G) - 6 m (18) 

g h (G) - G 6 (19) 

so that 

g r (G) - G m ' 6 (20) 

where n, m, and s are constants. 

EXTENSION TO ANISOTROPY 

The direction of grain growth or solidification at each point in a SC or 
DS solid can be characterized by a field of unit vectors d-j(x|<). The 
mechanical behavior at each point must then depend not only on the stress and 
deformation history at the point but also on the local preferential direction. 
This requires that dependence on d-j be included in F and G in equation 
(8). However, as the sense of d^ is immaterial, the dependence is properly 
taken in terms of the product d^dj. Thus, we replace equation (8) with 

F ( r ij’ d i d j> and G ( a ij» d i d j^ (21) 


As indicated in appendix B, the theory of tensorial invariants (refs. 8 
to 10) requires that, for form-invariance under arbitrary rigid-body rotations, 
F and G must be expressible in terms of the principal invariants of their 
respective tensorial arguments and invariants involving various products of 
these tensors. Here, as argued in appendix B, we take the functions F and 
G as depending on the subset of these invariants. 
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quality 


for F, and 


1 ^ 
! 1 * 7 x i J r ji 

l Z * d i d j E jk E ki ► 
’3 *7 d i d j I j1 



*7 a ij a j 1 

‘ d i d j a jk a ki 


f 



( 22 ) 


(23) 


(24) 

(25) 


As in the fully isotropic development we have sought generalizations of a 
von Mises type theory and have, therefore, restricted our choice of invariant 
expressions to those quadratic in stress. Analogous to the earlier develop- 
ment, K denotes the threshold (Bingham) shear stress transverse to the pref- 
erential material direction and K d denotes the same for shear along the 
material direction (fig. 1). For K * K d , indicating no difference in shear 
strength across and along the direction d^, equations (24) and (25) reduce 
to their isotropic counterparts (eqs. (11) and (12)). 

As before, the flow and evolutionary equations are obtained from equations 
(5) and (6), this time by making use of equations (7), (13), and (22 to 25). 
Again, the details are reserved for appendix A. The resulting flow law is 
given by 

• f < F * 


(26) 


£ 1 } M 72 - 1 VWkl * d k d 1 £ jk - 2 d iVki (s ij * d 


iv) 


and the evolutionary law by 
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*U * 9^57 - Rg r (G) 


a U + 


(I ' ') ("' 


ORIGINAL PAGE FS 
OF POOR QUALITY 

d k a ki + d k d 1 a jk 


-| d tVki (t ij + d i d j>) 


( 27 ) 


Note that e.^ « 0 and a. . « 0, the former indicating incompressibility of 
the inelastic deformation and the latter confirming the deviatoric nature of 
a-jj. As before, when K = equations (26) and (27) reduce to those of the 
isotropic case (eqs. (14) and (15)). 

Recall that the functions f, g and g^ and the material constants K, 
y, R, and H are determined just as in the isotropic case. Determination of 
K(j or alternately the ratio K/Kj is discussed in the following section. 


APPLICATIONS 


We first consider applications of the foregoing theory to the cases of 
homogeneously stressed elements in pure shear, transverse (fig. 1 (a)) and 
longitudinal (fig. 1(b)) to the preferential material direction. In each 
case, the coordinate direction is aligned with the material direction, 

i.e., £l = M. 0, 0 ). For creep in transverse shear (fig 1(a)), we have 


°11 


°22 = °33 = °12 * °13 = 0 

► 


^22 = x = const. J 


( 28 ) 


and 


“11 “ °22 " a 33 = a 12 " “13 = 


°23 = s 4 0 


> 

> 


From equation (24) we obtain for F 



( 29 ) 


( 30 ) 


which depends only on the parameter K. From equation (26), the flow law for 
the shear rate component ^3 = r tr is 


wr 


tr 


- s y 


-1 n 


- 1 
J 


(t - 5.) 


( 31 ) 
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where we have made use of the first of equations (17). 
longitudiral shear (fig. 1(b)), we write 

°11 ’ °22 * °33 - °13 - °23 * 0 


&12 ■ t ■ const. 


Now for creep under 




> 


j 


( 32 ) 


and 


a ll “ a 22 " a 33 “ “13 " a 23 “ 0 


a 12 “ s ^ ® 


(33) 


This time F from equation (24) is 


F 



(34) 


which now depends only on K .. The flow law for yi? 2 y , 0 likewise ob- 
tained from equation (26) and is ’ l 


UY 
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- ?) 


(35) 


With F » 0 and s»s *0 (corresponding to the initial stage of a creep 
test) we obtain by dividing equation (35) by equation (31) and solving for 


K 2 /K d 2 , 


7 

*d 



( l/n + l) 


(36) 


As the parameter n is assumed known, the ratio K/K^ is defined by equa- 
tion (36) in terms of the ratio of creep strain rates along and transverse to 
the preferential material direction. Simple shear tests of this type can, in 
principle, be used to determine the ratio K/Kj; however, a more practical 
method on the basis of uniaxial tests will be suggested in the following 
paragraphs. 

Next consider an application of the theory to that of uniaxial plane 
stress as depicted in figure 2. The plane stress element lies in the plane 
X3 « 0 with the tensile stress applied along xj. The preferential material 


v. ' 
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direction Is constant throughout the body and. In this case, wakes an angle f 
with the xj axis (fig. 2). The stress components are 

. ■n 

°2Z m °33 " °12 “ °13 " °23 “ 0 


with ORIGINAL PACE IS 
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■ o ^ 0 


(37) 


The unit vector d denoting the material direction Is 

d « (cos f, sin f , 0 ) 


(38) 


F and 6 from equations (24) and (25) become 

,2 


F . 

3IC 


1 ♦ 


- 3 ^ - 1 ^ ^4 cos 2 * ♦ sin 2 * - j (2 cos 2 * - s 1 n 2 f ) 2 j 


- 1 


(39) 


and 


3K‘ 


a 


^4 cos 2 * + sin 2 * - (2 cos 2 * - sin 2 *) 2 | 


(40) 


where a - is the uniaxial component of the internal stress The 

governing equations for the extensional strain rate c • cjj are given by equa- 
tions (26) and (27) as: 


1 


e ■ T" f ( F ) ( o - o) 


(4 cos 2 f ♦ 3 cos 2 * sin 2 * ♦ sin 2 f ) 


a ■ 


3H 


- Rg r (G), 


KH 

1 + ^^7 " ^ cos?f + 3 cos2f s1 ° 2f * 


sin 2 *) 


(41) 


(42) 


As expected, the in-plane shear strain rate y - Yj^ not generally zero and 
is given by: 
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y ” o) ^“ ^ 


The shear strain rate is zero when the preferential material direction Is 
0* or 9 - 90* end when the material Is Isotropic, l.e., K/K<j - 1. 

With f m 0* and f ( F ) - F n , equations (39) and (41) become 


3 vi F 0 (!^2‘ ) (o 



Under constant stress creep conditions with F >> 0, the ratio of initial 
creep rates corresponding to 9 * 0 * and 9 « 90* is given by 



2 2 

solving for K / we get 
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As the material parameter n is presumed known, equation (49) allows the 
determination of the ratio K/Kh from uniaxial tests with stress directed 
along and transverse to the grrin growth or solidification direction. Uniaxial 
tests conducted with the applied stress at arbitrary angles to the preferential 
material direction will provide information on wh^h assessments of the present 
theory can be made. 

For a complete elastlc-vlscoplastlc theory, needed for structural analy- 
sis, a compatible anisotropic elasticity formulation must be coupled with the 
present model. This will not be dealt with here but will constitute a topic 
of subsequent research. 


SUMMARY AND CONCLUSIONS I 

A constitutive theory has been presented for representing the anisotropic 
viscoplastic behavior of high-temperature alloys that have directional proper- 
ties resulting from controlled grain growth or solidification. The theory Is 
constructed by defining a vector field that Identifies the preferential direc- 
tion at each material point. This results In a locally transversely aniso- 
tropic model with allowance for spatially varying directional properties. The 
anisotropic theory Is based on the isotropic viscoplastic model of Robinson 
that has already been successfully applied In elevated temperature structural 
analysis. 

Application of the anisotropic theory is made to homogeneously stressed 
elements in pure shear with the shear direction taken transverse to and along 
the preferential material direction. These simple applications help to illus- 
trate the physical origin of the pertinent material parameters K and K^. 

Application is also made to a uniaxially stressed rectangular block in a state 
of plane stress with the spatially constant material direction making an arbi- 
trary angle with the stress direction. As expected, the results indicate that 
shear strain generally develops in the absence of shear stress, in other 
words, the principal axes of stress and strain are not in alignment as Is 
generally true under conditions of anisotropy. 

It is shown that the critical material parameter (or alternately the 
ratio K/K^) can be determined on the basis of uniaxial creep tests with the 
uniaxial stress direction along and transverse to the preferential material 
Jirection (grain growth or solidification direction). 

Tor a complete elasto-viscoplastic model, an appropriate elasticity form- 
ulation must be coupled with the present model. 
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APPENDIX A 

DERIVATION OF FLOW AND EVOLUTIONARY EQUATIONS 
Isotropic Case 


We first present the derivation of the flow and evolutionary laws for the 
fully isotropic case, i.e., leading to equations (14) and (15). From equation 
(5) we write: 


ao ao dF aJ 2 aS kt 

c . .• » 


ij a °ij TF dJg aS^ 3o,.j 


(1A) 


Making use of equations (7) to (12), we have 


and 


30 K x / r \ 

W f(F) 

(2A) 

dF 1 

(3A) 

3 0 2 


85 u ' 

(4A) 

3 ^kt 1 

"aoTT = 6 ki i i j “ '1 6 i j { ki 
1 0 

(5A) 


Substitution of equations (2A) to (5A) into equation (1A) leads directly to 
equation (14), 


2 -i0 


f(F)l, 


IJ 


(6A) 


Next, from equation (6) we write: 




-h 


an ^ dG 2 

aTd^aa kt 




where in addition to equations (2A) to (5A), we have 


an 

aE 


K 2 



9(G) 


(7A) 


( SA) 
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dG 



1 





-I 


ij 




and 


da 


kt 


1 


— ■- 6 * .6 , •• 6 . . 6 , 

ki ij 3 ij ki 


Combining these and using equation (13) gives equation (15), 

H 


*U-9?ET*1J- Rs r (S >*,J 


(9A) 

(10A) 

(11A) 

(1ZA) 

(13A) 


Anisotropic Case 


The derivation of equations (26) and (27) from equations (5) and (6) will 
now be outlined. From equation (5) we write: 


i j 


da da 



3F 8l l 

Tr T 85 ki 


. »F 3l 2 . »F 3l 3 

*^2 3$ ki sI 3 35 k» 




(14A) 


Using equation (7) together with equations (22) to (25) we have (without 
repeating terms already included in equations (2A) to ( 12A) ) 


dF 1 

d7~ = T 
31 1 k 


(15A) 


dl 


1 




(16A) 


dF 

tt 2 



(17A) 
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9^J * d jVki + d k d i I jk 


'fe'?) 


d 1 d j E ji 


3 *3 1 

^7 ■ 7 d ^ 


iJ ki 1 

9oT" * s ki 5 ij -7*1j*k» 

* J 


Substituting the appropriate terms from equations (2A) to (12A) and equations 
(15A) to (21A) into equation (14A) we get equation (26), 

2 v k.. - f (F) E -j j + ^ ^ ( d j d k r ki ' d k d i I jk " 2 d t d k £ ki* 6 ij + d iV) 


Finally, we write from equation (6) 

• aa l"aF 3l l . aF y -2 aF al 3 a fk t 

a ij aT TTY 3a kt TY aa k|[ al 3 aa Rt aa^ 


!!s& 


il 3 8d kt a^ 2 8a ki 9^3 8a kiJ 8a ij 


The terms in equation ( 2 3 A ) not previously stated are 


88 ij * - £ij 


fr - -fjVki * d k d i £ jk] 


3a ki 1 

da . ■ ki ij 3 ij ki 

* J 
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»G 1 

75; ‘ 7 




ta 


■ A 


u 


1J 


JLL /! 1 \ 

•VU'^7 


(2M) 

(28A) 

(29A) 


‘ d j d k 4 k1 ' Wjk 


'h v*S *7 


d 1 d j a j1 


3 -^3 1 

T7~ ’ 7 d j d < 


(30 A) 
(31A) 


(32A) 


Combining the appropriate terms in equation (23A) and again making use of equa- 
tion (13) we get equation (27), 


*ij • ? h TS7 ; ij ■ R, r (G) 


k c 


a ij * It?" 1 ) l d j d k a ki + d k d i a jk 


-? d i d k*ki (, u * d i d j>) 


(33A) 
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BASIS FOR SELECTION OF INVARIANTS 


It follows from the theory of algebraic invariants (refs. 8 to 10) that 
the scalar function F, specified in equation (21) as a function of the two 
symmetric tensors 




(IB) 


and 


d-ivj] 


(2B) 


is form-invariant under arbitrary rigid body rotations if expressed in terms of 
the invariants 


trD - trp 2 - trD 3 «= 1 


trS = 0, trS 2 , tr2 3 > 


tr D2 = trD 2 2, trD5 2 


J 


A set of nontrivial invariants extracted from equations (3B) is 

trS 2 , trS 3 , trDS and trOS 2 


(3B) 


(4B) 


Seeking generalizations of a von Mises type, we limit the dependence of F on 
combinations of the invariants in equations (4B) that are quadratic in stress, 
i .e. , 

’l trsZ ‘ 7 £ ij £ ji 


In particular, we take 


I 2 = trDS - 


(l trD5 J ■ (l d 1 d j £ jl)' 




(SB) 


(6B) 
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as expressed in equation (24). A parallel argument applies to the function G. 
According the equations (21) G depends on the symmetric tensors 

A - [a. . ] (78) 

ORIGINAL PAS~ ,3 
and ° F P00R WLITY 

O-C^dj] (88) 


and we are led, using arguments similar to the above, to express G in terms 
of the invariants 

h "jtrA 2 . -? a ij a ji "1 


$2 - trDA 2 = U 1 d j a J . k a ki 



(? trDA ) 2 * (? 


Specifically, we write 


> 


G = 




(9B) 


(10B) 


as given by equation (25). 
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figure 1. - ffomoqeneousl> stressed elements In pure sheer. 
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